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Abstract 

We study the transport properties of one-dimensional (ID) interacting 
Fermions through a barrier by numerically calculating the Kohn charge stiff- 
ness constant and the relative Drude weight. We find that the transport 
properties of the ID Hubbard model are quite different from those of the ID 
spinless Fermion model. For example, the presence of the attractive interac- 
tion between electrons in the ID Hubbard model actually suppresses the DC 
conductance, while a small repulsive interaction enhances the DC conduc- 
tance. These results show that the spin degree of freedom plays an important 

role in the transport properties of the ID interacting Fermion systems. 
PACS numbers: 71.30.-hh, 75.10.Jm, 73.20.Dx 
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Recently, there has been increasing interest in one-dimensional (ID) electron systems. 
On the one hand, the progress in ultrafine lithographic technology has made it possible 
to fabricate and experimentally study quantum wires in which only the lowest subband is 
occupied, approaching the true ID limit; on the other hand, there is a renewed interest in 
the Luttinger liquid because of the suggestion that the high temperature superconductors 
might be non-Fermi-liquid-like. It is well known that even a weak interaction in a ID 
Fermion system will drive the system to become a Luttinger liquid which exhibits quite 
different behavior from the Fermi liquid system. P,p!o| For example, Kane and Fisher have 



studied the transport properties in a ID spinless Luttinger liquid with a barrier (or a weak 
link). They find that an attractive interaction makes the system a perfect conductor while 
a repulsive interaction makes it an insulator. 

In this paper we present some results of the transport properties of ID interacting 
Fermion systems with a barrier. We study both the interacting spinless Fermion model 
and the ID Hubbard model with a barrier by calculating the Kohn charge stiffness constant 
using the Lanczos algorithm. In particular, we investigate the effect of spin on the transport 
properties motivated by an interesting statement of Anderson who argued that the theory 
of the transport in ID spinless Fermion may be irrelevant for a real electron system which 



has both charge and spin degrees of freedom. |Tl|] We find that our result for the spinless 
Fermion model agrees with that of Kane and Fisher qualitatively. But our result of the 
ID Hubbard model with a barrier is exact opposite of that of the spinless Fermion model. 
We find that in the case of the Hubbard model, the presence of the attractive interactions 
between electrons actually decreases the Kohn charge stiffness constant D^ as well as the 
relative Drude weight, in contrast to the ID spinless Fermion case where the presence of 
attractive interactions increases Dc and the relative Drude weight as long as the attractive 
interaction is not large enough to cause a phase transition to a phase-separated insulating 
state. In the case of repulsive interaction, we find that a small repulsive interaction in the 
ID Hubbard model with a weak link actually increases D^ and the relative Drude weight, 
again in contrast to the ID spinless Fermion model where the repulsive interaction decreases 



We first study an interacting spinless Fermion model, so called Heisenberg-Ising model, 
on a ID ring with a weak link 

H = -^ti(cjci+i + h.c.) + V^^n^ni+i. (1) 

i i 

Here ti = t for i ^ I and ti = t' representing the barrier (or weak link) site, and V is the 
nearest neighbor interaction. We assume the periodic boundary condition and thread the 
ID ring with a flux $ which is represented by a vector potential A = ($/L)ex, where L is 
the number of sites. We calculate numerically the ground state energy of the system E{^) 
as a function of flux, and consequently the Kohn charge stiffness constant 

Dc is a measure of the electronic conductivity of the system, and it vanishes if the system 
is an insulator. Kohn first suggested that the charge stiffness constant can be used as a 
quantitative measure of the Mott metal-insulator transition. [|T^-|TB 



Physically, the charge stiffness constant Dc is just the Drude weight in the real part of 
the optical conductivity ai{uj) in the long wavelength limit, 

ai{uj) = 2Tie^D,6{uo) + a\'\uj), (3) 



here we have set U = 1. Using second order perturbation theory, one finds, |16 



d. = ^<-t>-1e'^°'^'^'"^^^ (4) 

where < —T > is the expectation value of the kinetic energy in the ground state. ai{u) 
satisfies the well-known f-sum rule, 

ai{uj)duj = — < -T > . (5) 

-oo 1j 

We use the Lanczos algorithm to calculate the ground state energy. For a 14-site ring 
with 7 particles, the dimension of the Hilbert space is about 3400 and the Hamiltonian 



matrix has about 30,000 non-zero elements which can be handled easily on a workstation. 
We first consider a ID ring with the lattice translational symmetry, i.e., t' = t, in which most 
of the properties are known. This well-studied case provides us with a test example which 
we can check our computer program against and get some feelings about the finite size effect, 
etc. In Fig. la, we show the calculated ground state energy E{^) as a function of flux $ for 
a 14-site ring with 7 and 6 particles, respectively. Obviously, £'($) is a periodic function of $ 
with a period $o = hc/e. For the 7 particle system, the energy minimum is located at $ = 
as expected. This is true as long as the number of particles is an odd integer. But when there 
are even number of particles on a ring, such as in the 6-particle case, the minimum of E{^) 
is located at $ = $o/2, signaling that the system prefers a spontaneous persistent current 



which would generate a flux of $o/2- This is the well-known even-odd finite size effect. |jT7 
Fig. lb shows the calculated Kohn charge stiffness constant Df. as a function of the nearest 
neighbor interaction V with L=14, N=7, and t' = t. One can see that on the repulsive 
interaction side, approaches zero quickly as the repulsive interaction increases, signaling 
a transition to the Mott-insulator in this half-filling case. On the attractive interaction side, 
Dc also start to decrease, even faster than the repulsive case, at around V = —2t. This 
is the transition to a phase-separated insulator state. On the repulsive interaction side, it 
is well known that Dc and the effective mass have a discontinuous jump at = 2t in the 
thermodynamic limit. [|l^ As \^ ^ 2t — 0, approaches a nonzero value 1/4, while = 
as V ^ 2t + 0. Our finite size calculation yields ~ 0.285 at = 2t. At V=0, Dc in the 
thermodynamic limit is l/vr ~ 0.3183. Our calculation gives — 0.3210. The finite size 
calculation usually overestimates Dc because of the periodic boundary condition. 

In Fig. 2a, we show Dc as a function of V for a 5-particle 16-site system with t' =t, 0.75t, 
and 0.25t, respectively. Fig. 2b shows the relative Drude weight (RDW) which is defined 
as the ratio of the Drude weight Dc and the total spectral weight Afot = ai{uj)duj/2TTe'^ 
for the same system. The first thing we can see from Fig. 2 is that the presence of a weak 
link does not seem to affect the phase-separation transition much. The weak link actually 
makes the phase-separation-transition sharper in a finite size system. The transition is 



depicted better in the relative Drude weight plot (cf. Fig. 2b). It is obvious that any 
theory which assumes a uniform phase such as that of Kane and Fisher ^ is only applicable 
above the threshold Vc = —2t. Below Vc, Dc is zero and the system is an insulator in 
the thermodynamic limit. The second thing we observe from Fig. 2 is that the attractive 
interaction does enhance the conductivity, especially when t'/t is small, that is quite different 
from the result of ID Hubbard model which we'll discuss later, while the repulsive interaction 
suppress the conductivity. This is in qualitative agreement with the analytical studies of 
Kane and Fisher [0], and Apel and Rice. However, due to the finite size effect and 
the phase-separation transition at = —2t, presently we are not able to confirm Kane 
and Fisher's specific prediction that any repulsive interaction will drive the system to 
insulating while any attractive interaction will make the system a perfect conductor. We find 
that Dc and the relative Drude weight remains finite even for very large positive V in a finite 
system. The third point we learn from Fig. 2 is that when t' = t, the Drude weight almost 
exhausts all the spectral weight Afot as long as V is above the phase-separation threshold 
Vc = —2t. The relative Drude weight drops quite rapidly once V is decreased below V^. 
When t' < t, the RDW is usually smaller than the t' = t case (cf. Fig. 2b). This is due 
to the pinning of the charge density wave by the weak link, which not only reduces Dc, but 
also reduces the RDW. 

Next we consider the Hubbard model on a ID ring with a weak link, 

icr i 

Here again ti = t' and ti = t for i ^ 1, and o" =t, | is the spin index. In Fig. 3 we show the 
calculated Dc and the relative Drude weight as functions of the interaction U for t' =t, 0.5t, 
and O.lt, respectively. We find that the transport properties of the ID hubbard model are 
quite different from those of the spinless Fermion model. In Hubbard model, the presence 
of the attractive interaction between electrons actually decreases Dc and RDW (cf. Fig. 
3), in contrast to the spinless Fermion model (cf. Fig. 2 and Refs. |]7|,pr8[|). The result of 
the repulsive interaction is equally surprising. The presence of a small repulsive interaction 
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{U > 0) increases Dc and RDW. So the maximum of Dc actually shifts from U=0 for t' = t 
to some small positive UQ{t') for t' < t (cf. Fig. 3). f/o(^') increases as t' decreases. For 
t' = O.lt, Uq ~ t. Since the only major difference between the Hubbard model and the 
spinless Fermion model is the spin, it is reasonable to suspect that the spin does play a 
very important role in the transport of ID interacting Fermions, which is in agreement with 
Anderson's conclusion, fill 



The fact that in Hubbard model, Dc decreases as one increases the strength of the at- 
tractive interaction can be understood from the following physical picture. In the attractive 
interaction region, electrons form pairs which can be treated as hard-core bosons, |jl9|, and 



these bosons can hop via virtual ionization. A straight forward second-order perturbation 
calculation yields that the pair hopping amplitude is 2t^ /\U\ and the nearest neighbor re- 
pulsion between the bosons is also 2t^/\U\ in the t' = t case. Physically, this is easy to 
understand because the virtual ionization to a nearby site lowers the energy, but it is only 



possible if the nearby site is empty. ||20| 



Comparing our results with that of Ref. 7 and 8, we find that although our results in the 
spinless case is in qualitative agreement with their results, there is significant discrepancy 
between our results and that of Ref. 7 and 8 in the spinful Fermion case. The Hubbard 
model we studied has SU(2) symmetry, so is fixed at 2 (we use the same notation g„ as 
in Ref. 7). This means, according to Ref. 7 and 8, that the spin part should behave like the 
noninteracting case while the charge part should act like the spinless Fermion. This is in 
disagreement with our finding that the presence of a small repulsive interaction in Hubbard 
model with a week link actually increases DC conductance which is quite different from the 
results of spinless Fermions. One possible explanation of our results is that the presence 
of repulsive interaction increases the antiferromagnetic correlation which may enhances the 
transport across the weak link. Some of these physics is apparently missing in Ref. 7 and 
8. For example, Ref. 8 starts from a model which only has interaction g2 and completely 
neglects the 2kp component gi. Here gi and g2 are standard "g-ology" notations as in Ref. 
9. 
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Finally we mention that although all of the results shown above are for t' <t, we have 
also studied the case of t' > t. In this case, we find that Dc decreases as t' is increased from 
t, in agreement with the picture that the t' ^ t link (either t' < t or t' > t) serves as a 
pinning center of the charge density wave in the repulsive interaction region. 

In conclusion, we have studied the transport properties of both the ID spinless Fermion 
model and the Hubbard model with a weak link by calculating the Kohn charge stiffness con- 
stant Dc and the relative Drude weight using Lanczos algorithm. We find that the presence 
of the attractive interaction between electrons in ID Hubbard model actually suppresses the 
DC conductivity, which is in contrast to the ID spinless Fermion model where the attractive 
interaction enhances the DC conductivity as long as the attractive is not large enough to 
cause a phase transition to a phase-separated insulating state. In the repulsive interaction 
region, we find that a small repulsive interaction in ID Hubbard model actually increases 
Dc and the relative Drude weight, again in contrast to the ID spinless Fermion model where 
the repulsive interaction decreases D^. These results show that the spin degree of freedom 
plays an important role in the transport properties of the ID interacting Fermion systems. 
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FIGURES 

FIG. 1. (a). Shows the calculated ground state energy E{^) of the Heisenberg-Ising model as 

a function of the flux $ for N=7 (solid) and 6 (dashed), respectively. Here L=14, V=2t, and t' = t. 
The unit of E is t. (b). Shows the calculated Dc as a function of interaction V for the N=7 system. 

FIG. 2. Shows the calculated (a) Dc and (b) the relative Drude weight as a function of inter- 
action V of a 5-spinless-Fermion system on a 16-site ID ring with t'=t (solid), 0.75t (dashed), and 
0.25t (dot-dashed), respectively. 

FIG. 3. Shows the calculated (a) Dc and (b) the relative Drude weight of a ID Hubbard model 
as a function of interaction U of a 6-Fermion system on a 8-site ring with t'=t (solid), 0.5t (dashed), 
and O.lt (dot), respectively. The total spin Sz = 0. 
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